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Summary  of  overall  progress 


The  analysis  of  the  advantages  and  disadvantages  of  the  homomorphic 
filtering  technique  applied  to  solve  the  exact  inverse  scattering  problem  has 
continued.  There  have  been  problems  with  strongly  scattering  objects  due  to 
phase  unwrapping.  Work  has  continued  to  more  carefully  compare  the 
reconstruction  of  V'F  from  the  backpropagation  of  simulated  data  with  that 

predicted  for  the  direct  problem,  from  the  same  data. 

We  note  that  piiase  unwrapping  is  avoided  with  the  differential  cepstrum; 
we  believe  that  there  might  be  an  interesting  connection  between  this  and  a 
novel  differential  formalism  for  the  Rytov  approximation. 

This  report  focuses  mainly  on  new  developments  in  modelling  scattering 
from  nonlinear  media.  In  anticipation  of  the  next  phase  of  this  project,  namely 
to  consider  the  inverse  scattering  problem  for  nonlinear  medium,  we  have  been 
studying  the  writing  of  structures  into  nonlinear  media,  specifically 
photorefractive  crystals.  A  paper  on  this  has  been  rewritten  following 
submission  to  JOSA-B  and  has  now  been  submitted  to  Opt.  Comp.  &  Sig.  Proc. 
Both  experimental  and  theoretical  work  continues  toward  writing  structures 
with  intense  beams  and  reading  with  weak  beams,  the  partial  erasure  on  reading 
being  a  problem  without  fixing.  Methods  for  fixing  are  being  independently 
developed  within  the  group.  '  »  '  - 
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Activities 


During  the  last  four-month  period 

Next  July,  M.A.  Fiddy  is  chair,  and  F.C.  Lin  co-chair  of  the  Conference  on 
Inverse  Problems  in  Scattering  and  Imaging,  at  SPIE’s  1992  International 
Symposium  on  optical  Applied  Science  and  Engineering,  in  San  Diego.  There 
are  50  papers  from  aU  over  he  world  being  presented. 


Papers  and  conference  presentations 


Lin,  F.C.  and  M.A.  Fiddy,  "On  the  issue  of  the  Bom-Rytov  controversy:  I 
Comparing  analytical  and  approximate  expressions  for  the  one-dimensional 
case",  accepted  J.O.S.A.  A. 

Lin,  F.C.  and  M.A.  Fidd  "Optimization  of  the  self-pumped  phase  conjugation 
in  BaTiOj  for  optical  image  storage  and  readout",  submitted  to  Opt.  Comp,  and 
Sig.  Proc. 


There  are  also  2  papers  on  inverse  problems  to  be  presented  at  the  San  Diego 
meeting  in  July. 


Lin,  F.C.,  R.  McGahan,  A.  Alavi  and  M.A.  Fiddy,  "Quantitative  image 
recovery  and  restoration  from  scattered  field  data  at  lOGHz,"  to  be  presented  in 
the  Conference  on  Inverse  Problems  in  Scattering  and  Imaging,  at  SPEE's  1992 
International  Symposium  on  optical  Applied  Science  and  Engineering,  in  San 
Diego. 

McGahan,  R.,  F.C.  Lin  and  M.A.  Fiddy,  "Cepstial  filtering  for  recovery  of 
object  from  scattered  field  data,"  to  be  presented  in  the  Conference  on  Inverse 
Problems  in  Scattering  and  Imaging,  at  SPIE's  1992  International  Symposium 
on  optical  Applied  Science  and  Engineering,  in  San  Diego. 
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Progress  since  last  report 

1 .  Propagation  of  Light  in  One-dimensional  Half-space  Nonlinear  Media 

Because  of  their  large  nonlinearities  observable  with  the  usage  of  lower- 
power  cw  lasers,  artificial  Kerr  media  such  as  aqueous  suspension  of 
submicrometer  dielectric  spheres  may  be  found  useful  as  nonreciprocal  media 
for  switching,  image  amplification  or  as  storage  media  for  holographic  imaging 
[1].  We  are  interested  in  their  use  as  controllable  scattering  media  for  design 
and  synthesis  of  optical  components,  in  which  strong  scattering  phenomena  play 
a  role. 

Consider  a  laser  beam  incident  from  free-space  (region  0  for  z  <  0)  into 
the  one-dimensional  (1-D)  half-space  nonlinear  medium  (region  1  for  z  >  0) 
[Figure  1],  Suppose  that  the  incident  laser  beam  is  a  time-harmonic  plane  wave 
^i(koz  G)t)  radial  frequency  co  and  =  CD^/p^  where  Po  and  Eq  are  the 

free-space  permeability  and  permittivity,  respectively.  Also  assume  that  the 
nonlinear  medium  can  be  characterized  by  the  free-space  permeability  Po  and  a 
permittivity  ei(z)  =  ei[no(z)  +  An(z)]^.  The  spatial  distribution  of 
inhomogeneities  in  the  medium  modulates  locally  the  linear  refractive  index 
no(z)  [1]  while  the  effects  of  self-actions  of  light  in  the  medium  gives  rise  to  the 
nonlinear  term  An(z)  [2].  For  example,  in  artificial  Kerr  media  such  as  liquid 
suspensions  of  submicrometer  dielectric  particles,  the  electrostrictive  effect 
yields  the  intensity-dependent  contribution  to  the  nonlinear  term,  i.e., 

An(z)  =  |^l4',(z)P  (1) 
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Figure  1. 


where  =  y  eJzq  and  n2  and  M^i(z)  are  the  nonlinear  (Kerr)  constant  and  the 

total  field  in  region  1  (the  nonlinear  medium),  respectively.  For  liquid 
suspension  of  submicrometer  polystyrene  latex  spheres  with  0.234  pm  in 
diameter,  the  typical  value  of  n2  is  about  10'^  cm^AV  [3].  Therefore,  for  the 
moderate  laser  intensity  from  1  to  100  W/cm^,  An(z)  can  be  taken  as  a  small 
perturbation  to  no(z). 

In  the  1-D  scattering  theory,  the  total  fields,  To(z)  and  'Fi(z),  in  rep' on  0 
and  1  satisfy  the  scalar  wave  equations 

i-^^+kJ%(z)  =  0;  z<0  (2) 

dz 

and 


4 


d^4^,(z)  ,2 


dz 


f^+kiH',(z)  =  -k; 


*^0(2)  +  r;fl4"i(z)l 


2n, 


-  1  H'lCz);  z>0  (3) 


where  kj  =  .  For  brevity,  the  time-harmonic  factor  e  is  omitted 

from  the  total  fields.  The  term  on  the  right-hand  side  of  Eq.  (3)  is  the  so-called 
secondary  source. 


1.1.  The  Bom  Method  for  the  1-D  Half-space  Nonlinear  Case  with  no(z)  =  1 
and  El  =  Co  =  1) 

1  ik  Iz  -  z^l 

In  this  case,  the  free-space  Green's  function  G(z,  z')  =  ® 

should  be  used  in  the  inhomogeneous  Fredholm  equations  of  first  kind  for  To(^) 
and  'Fi(z)  [4],  namely, 


j 


z>0  (5) 


0 

where  Tq  (z)  and  4^,  {z),  which  are  also  equal  to  the  incident  plane  wave 
for  this  case,  are  the  unperturbed  fields  in  regions  0  and  1,  respectively.  Since 


the  contribution  due  to  -k3  0214^1  (z)r  +  l4'i(z)l 


n| 


4^1  (z)  to  the  scattered  fields 


is  very  small,  we  can  apply  a  perturbation  method  to  expand  the  total  fields  of 
Eqs.  (4)  and  (5)  in  terms  of  the  unperturbed  fields,  namely, 


5 


^^(z)  =  tIV)  -  kj  Jdz'G(z,z')[njlH'nz')  +  H-'.'tz')  +  +  ...I" 

0 

+  +  'Flstz')  +  'J'utz')  +  ...1'*]  ['FiTz')  +  'f\^'^z')  +  T^i^s1:z')  +  ...] 


m  =  0,  1  (6) 

where 'F  is  tz)  (n  =  1,  2,  ...)  is  the  nth-order  Bom  approximated  scattered  field 
in  region  1.  After  lengthy  calculation,  we  obtain  the  first-order  and  the  second- 
order  Bom  approximated  scattered  fields  in  region  0  and  1 : 


and 


z  <  0, 

(7a) 

z  <  0, 

(7b) 

'f1sV)  =  (-  T+  i  Tkozje‘^®^ 

z  >  0, 

(8a) 

z  >  0. 

(8b) 

1.2.  Exact  Solutions  of  HcCz)  and  Ti(z)  for  the  1-D  Half-space  Nonlinear  Case 
with  no(z)  =  1  and  Ei  =  Zq  (nj  =  1) 

The  closed-form  solutions  for  To(z)  and  T'i(z)  can  be  written  as 

'Fo(z)  =  +  Rs  z  <  0  (9) 

and 

YiCz)  =  T,  e'*'’";  z>0.  (10) 

The  real  reflection  and  transmission  coefficients,  Rj  and  T^,  can  be  derived 
from  the  boundary  conditions: 
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%{0)  =  T.(0) 

and 

d%(z)  _  d%{z) 

dz  z  =  0  ^  z  =  0. 

These  boundary  conditions  can  be  applied  for  the  acoustic  case  or  the  TE  case  in 
the  electromagnetic  wave  theory.  Since  the  wavenumber  in  the  nonlinear 
medium  is  boundary  conditions  yield  two  equations  for 

Rj  and  Tj,  i.e., 

l  +  Rs  =  T,  (13a) 

and 

1 -R,  =  (l  +  (1 3b) 

from  which  we  obtain 

njT’  +  4T,-4  =  0.  (14) 

Eq.  (14)  has  one  real  root  given  as 


and  two  complex  roots  which  are  complex  conjugate  each  other.  The  real  root 
of  Eq.  (14)  gives  the  correct  approximation  for  the  total  fields  when  n2  «  1, 

namely. 


^2  3n2 


and 


(16) 
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nj  .n2 


3n?  Sn? 


„  2 

"2  ,  2 


4^1  (z)  ~  1  -  +  i  X  koz  +  -  i  koz  -  ^  ^oZ 


...le''^.  (17) 


Note  that  Eqs.  (16)  and  (17)  are  exactly  equal  to  the  Bom  approximated  fields 
derived  from  Eqs.  (4)  and  (5). 

1.3.  The  Bom  Method  for  the  1-D  Half-space  Nonlinear  Case  with  no(z)  =  1 
and  £i  >  Eo 

In  this  case,  the  Green's  functions,  Goo(z,  z^,  Gio(z,  z'),  Goi(z,  z'),  and 
Gii(z,  z')  which  are  governed  by  the  equations 


d  u,v,(z,  z') 

- ^^2 - -  +  Is  Gp<,(z.  z)  -  6„5(z .  zO; 


(p,q  =  0,  1)  (18) 


are  used  in  the  inhomogeneous  Fredholm  equations  of  first  kind  for  4o(z)  and 
Ti(z).  Note  that  5pq  is  the  Kronecker  delta  and  5(z  -  z')  is  the  delta  function. 
Physically,  Gp^^(z,  z")  describes  how  the  field,  generated  by  a  point  source  at  z' 
(in  region  q),  behaves  at  the  observation  location  z  (in  region  p).  From  the 
boundary  conditions  at  z  =  0: 

G,J0,  z')  =  G,p(0,  zO  (19) 


dGpp(z,  z') 
dz 


z  =  0 


dGpp(z,  z') 

^  z  =  0 


where  p,  q  =  0,  1  and  q  p,  we  obtain 


-iko(z  +  z'Y 


z&z'<0,  (21) 
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G,„(z,  z-)  =  -^-T(k,)e‘'‘‘V''^"'; 

z'<0<z,  (22) 

G„(z,  z-)  =  i^[e‘''''^-^''-R(k.)e““'^^^'>]; 

z&z'>0,  (23) 

and 

G„,(z,  z')  =  T(k,)e'“‘“"e“‘'"': 

z  <  0  <  z'.  (24) 

Note  that  the  coefficients 

Icq  “  k] 

(25) 

and 

2ko 

(26) 

are  the  same  as  the  reflection  and  the  transmission  coefficents  in  the  three- 
dimensional  case  for  a  TE  wave  normally  incident  on  a  dielectric  half-space 
with  permittivity  Ej.  Owing  to  the  symmetry  of  the  scattering  geometry  in  the 

dielectric  half-space  case,  the  Green's  functions  satisfy  the  symmetry  properties 
[5],  namely, 

Gpq(z,  z)  =  G^(z,'z);  (p,  q  =  0,  1).  (27) 

Therefore,  the  total  field  T'p(z)  in  region  p  (0  or  1)  can  be  expressed  as 


T„{3)  =  t'7z)  -  k 


dz'Gpi(z,z) 


n 


n\ 


4nt 


(28) 


where  in  this  case,  the  unperturbed  fields,  T'o°\z)  and  4^i^°\z),  are  equal  to 
1  Ic  z 

and  T(kj)e  '  ,  respectively.  After  applying  the  perturbation  method  to  expand 
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Eq.  (28)  into  the  Bom  series  simliar  to  Eq.  (6)  and  using  Eqs.  (23)  and  (24),  we 
obtain  the  first-  and  the  second-order  Bom  approximated  scattered  fields 


IBA 


-ikoz 


2BA  ^^7  7  -ikrtZ 

'J'os  (z)  =  -irT’(k,)e 


'f'ls  (z)  =  [-  TT^(ki)  +  i  j;^T\k,)k,z]e‘'‘‘^ 
and 


2BA 

'i'is(z)  = 


3n2  7  3n2  x  <■  -» 

-I5-T'(k,)  -  i  S5f-T«(k,)k,z  -  j;|T^k.)k? 


ikfZ 

h  A  • 


z  <  0,  (29a) 

z<0,  (29b) 
z  ^  0,  (30a) 


z^O.  (30b) 


1.4.  Exact  Solutions  of  %(z)  and  'Fi(z)  for  the  1-D  Half-space  Nonlinear  Case 
with  no(z)  =  1  and  Si  > 

Similar  to  Eqs.  (9)  and  (10),  Ho(z)  and  'F2(z)  can  be  written  as 

%(z)  =  +  Rj  z  <  0  (31) 

and 

T.cz)  =  Tj z>o  (32) 

where  k^  =  kj  -f-  Using  the  same  boundary  conditions  as  in  Eqs. 

(11)  and  (12),  we  can  derive  two  relations  for  and  from  w-hich  we  obtain 

"r  Td’ +  Td  -  4  =  0.  (33) 

Eq.  (33)  has  also  one  real  root  given  as 


10 


and  two  complex  roots.  Thus,  with  Eq.  (34)  and  the  condition  that  n2  «  1,  the 
total  fields  can  be  approximated  as 


Eqs.  (35)  and  (36)  are  exactly  equal  to  the  Bom  approximated  fields  in  Eqs. 
(29a)  to  (30b). 


2.  Image  Reconstructions  for  RADC  Bistatic  Scattering  Measurements 

In  a  previous  report  (October  1,  1990  to  April  8,  1991),  we  applied  the 
far-field  filtered  backpropagation  algorithm  to  reconstruct,  from  the  far-field 
scattered  data  taken  at  the  RADC  bistatic  scattering  experiments,  the  images  of 
two  empty  annular  cardboard  cylinders  with  8.5  cm  and  15  cm  in  diameter  and 
0.4  cm  in  thickness.  These  reconstructions  are  accomplished  by  backpropagating 
36  sets  of  10  GHz  scattered  field  data  'BfCp*,  kj^^)  into  the  scatterer  domain. 

Each  set  of  data,  containing  181  data  points,  corresponds  to  one  viewing 
direction;  the  viewing  angle  ranges  from  0°  to  350®  with  10°  increments.  Each 
data  point  within  each  set  of  data  corresponds  to  one  scattering  angle;  the 
scattering  angle  ranges  from  0°  to  180°  with  1°  increment.  Before  using  the  far- 
field  backpropagation  algorithm  resulted  from  Eq.  (23)  of  the  previous  report, 
the  scattered  data  is  interpolated  into  the  cartesian  grid  on  the  k-domain  in 
order  to  apply  the  fast  Fourier  transform.  The  interpolation  module  (the  IMSL 
subroutine  "IQHSCV")  involves  a  smoothing  process.  Also  the  summation  over 
all  viewing  angles  (using  Eq.  (23)  of  the  April  1991  report)  combines  both  the 
along-  and  the  across-viewing-direction  structures  of  V(p)T^(p)  for  all  views 
(where  V(p)  is  the  scattering  function  and  T(p)  is  the  total  field  in  the 

scatterer).  Therefore,  both  the  along-  and  the  across-viewing-direction 
structures  of  V(p)T'(p,  ko^^)  for  each  view  with  the  viewing  direction  have 

been  artificially  symmetrized.  Although  the  far-field  backpropagation  algorithm 
provides  a  way  of  fast  reconstruction,  the  brute-force  method  is  necessary  to 
integrate  Eq.  (20)  of  the  April  1991  report  in  order  to  reconstruct  the  true 
along-  and  the  across-viewing-direction  structures  of  V(p)T(p,  ko0o)- 

Suppose  that  the  scattered  field  data,  ^^(pj,  ko^o)’  w'ere  measured  in  the 
far-field  at  positions  p^  =  p,  (5  coscp*  +  ^  sintpj  after  the  scatterer  had  been 


illuminated  by  a  time-harmonic  plane  wave  of  wavenumber  propagating  in 
the  direction  ^  cos(p„  4-  ^  sincp^,.  In  2D  diffraction  tomography,  we  have 
shown  [6]  that  the  scattered  field  ^s(p*,  is  written  as, 


ij2 

koOo)  =  X  Jd'p'  H'„'’(k„lp.  -  p-'l)V(p'mp',  k„^„)  (37) 

Q 

where  Hq  ^(kolp,  -  p  'I)  is  the  zeroth  order  of  the  Hankel  function  of  first  kind 
and  the  integration  is  over  the  cross  section  Q  of  the  scatterer.  Under  the  far- 
field  approximation,  we  have 

Hj'lkolp,  -  ri)  =  K/tT-  '  (38) 

and,  thus 

'i'.(p..k<,P.) »  y  J'd’p'e-“k.f'.-rv(p')'i'(p',k„^,).  (39) 

By  defining  a  function 

f(^  k„&j  =  V(p)'i'(p,  koP„)  (40) 

then,  Eq.  (39)  can  be  rewritten  as 


K  =  ^kQ(cos(p,  -  coscpo)  +  $)k^(sin(p,  -  sincpo).  (42b) 

After  taking  the  inverse  Fourier  transform  of  Eq.  (41),  we  obtain 

OQ  OO 

f(p.a.)=  JdK.  JdK,e‘(K.x  +  K,y)p^^(K,_K,).  (43) 

-OO  -OO 

In  fact,  the  limits  of  integration  in  Eq.  (43)  can  only  be  from  to  2ko 
because  the  measured  data  F(p^(K3j^  Ky)  are  confmed  within  a  disc  centered  at 

the  origin  and  having  a  radius  2ko  on  the  k-domain.  If  we  transform  the 
integration  dK^dKy  to  dtp^dtp,  by  using  the  relations: 


and 


=  lq)(cos(p,  -  costpo) 
Ky  =  ko(sin(p,  -  sin(p„). 


(44a) 

(44b) 


Eq,  (43)  becomes 

2  2k  2k 
d<p„Jd(p. 

which  is  the  same  as  Eq.  (20)  of  the  previous  report  that  is  used  to  formulate 
the  far-field  filtered  backpropagation  algorithm.  Note  that  the  dependence  of 
the  viewing  angle  in  4^(p,  k^^o)  is  averaged  out  in  Eq.  (45).  However,  if  only 

one  viewing  angle  is  considered  in  Eq.  (45),  we  can  reconstruct  the  along-  and 
across-viewing-direction  structures  of  V(p)4^(p,  ko^o)- 


V 1  -  F,pJk<,(^.-&.)).  (45) 


V(p)'{'(p)  =  ^  f 


In  Figures  2  to  7,  the  2-D  grey-level  images  are  the  modulus  plots 
reconstructed,  from  two  sets  of  simulated  and  two  sets  of  measured  scattered 
field  data,  over  a  spatial  area  64  x  64  cm^  in  the  object  domain  while  the  1-D 

linear  plots  are  real  parts  of  the  corresponding  images  either  along  or  across  the 
viewing  direction.  For  comparison,  the  actual  size  of  each  cylinder  is  indicated 
by  the  tip  of  the  arrowhead  on  each  2-D  grey-level  plot.  Figures  2a  and  3a  are 
the  reconstructions,  from  simulated  scattered  field  data,  of  cross  sections  for 
solid  cylinders  both  having  V  =  0.03  and  with  4.25  cm  and  15  cm  in  radius, 
respectively.  Figures  2b  and  3b  (2c  and  3c)  are  the  along-viewing-direction  (the 
across-viewing-direction)  plots  of  the  real  parts  for  the  objects  corresponding  to 
Figure  2a  and  3a  along  the  line  y  =  0  (x  =  0).  Figures  4a  and  5  are  the 
reconstructions,  respectively  from  measured  scattered  field  data  for  incident 
angle  (or  viewing  angle)  at  0°  and  240°,  of  cross  section  for  the  annular 
cardboard  cylinder  having  a  radius  4.25  cm,  thickness  0.4  cm,  and  filled  with 
styrofoam  of  V  =  0.03.  Figure  4b  and  4c  are  the  along-  and  the  across-viewing- 
direction  plots  of  the  real  parts  for  the  object  corresponding  to  Figure  4a  along 
the  lines  y  =  0  and  x  =  0,  respectively.  Figures  6a  and  7a  are  the 
reconstructions,  respectively  from  measured  scattered  field  data  for  incident 
angle  at  0°  and  90°,  of  cross  section  for  the  empty  annular  cardboard  cylinder 
having  a  radius  15  cm  and  tliickness  0.4  cm.  Figures  6b  and  7b  (6c  and  7c)  are 
the  along-viewing-direction  (the  across-viewing-direction)  plots  of  the  real 
parts  for  the  objects  corresponding  to  Figure  6a  and  7a  along  the  line  y  =  0  (x  = 
0).  Note  that  after  the  filtering  algorithm  has  been  applied  to  remove  noise 
from  experimental  data,  we  are  able  to  obtain  consistent  reconstructions  for  two 
different  viewing  directions;  comparing  Figures  6b  and  6c  with  Figures  7b  and 
7c. 
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Figure  6b. 
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Figure  7a. 


Figure  7b. 
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Figure  7c. 
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Summary  of  work  directions  in  the  next  period 


The  following  areas  will  be  pursued  in  the  next  few  months: 

i)  Homomorphic  filtering  techniques  applied  to  the  inverse  scattering 
problem,  for  simulated  and  real  data. 

ii)  Inversion  of  data  for  the  recovery  of  nonlinear  permittivity  profiles, 
with  the  expectation  that  we  can  design  permittivity  structures  in  nonlinear 
media  such  as  photorefractives. 
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